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Abstract. Results of a ground-based optical monitoring 
campaign on NGC 5548 in June 1998 are presented. The 
broad-band fluxes {U,B,V), and the spectrophotometric 
optical continuum flux Fa (5100 A) monotonically decreased 
in flux while the broad-band R and / fluxes and the inte- 
grated emission-line fluxes of Ha and H/3 remained con- 
stant to within 5% . On June 22, a short continuum flare 
was detected in the broad band fluxes. It had an ampli- 
tude of about ^18% and it lasted only ^90 min. The broad 
band fluxes and the optical continuum F\(5100 A) appear 
to vary simultaneously with the EUV variations. No reli- 
able delay was detected for the broad optical emission lines 
in response to the EUVE variations. Narrow H/3 emis- 
sion features predicted as a signature of an accretion disk 
were not detected during this campaign. However, there is 
marginal evidence for a faint feature at A ~ 4962 A with 
FWHM~6A redshifted by Av~1100 kms'^ with respect 
to H/3, 



able on timescales of a few days until years. Variations on 
timescales of hours and even less have been observed in 



-'narrow • 
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1. Introduction 

The ultraviolet and optical continuum and the broad emis- 
sion line flux of Seyfert 1 galaxies are known to be vari- 
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X-rays (e.g. Barr fc Mush otzky |1986| ; Leighly et al. |1996 
George et al. |l9984|l998b| ) 



To study the physical processes which are responsible for 
the observed spectral energy distribution of an active galac- 
tic nuclei (AGN) multiwavelength monitoring campaigns 
have proven to be an excellent tool (cf. Peterson 1993; Net- 
zer & Peter son |1997 for a review; Marshall et al. 1997; 
Nandra et al. |l998[ Edelson et al. |2000| ). Thus, over the last 
decade, several large space-based and ground-based mon- 
itoring programs have been undertaken for nearby AGN, 



such as NGC 5548 (Clavel et al. 1991 Peterson et al 



1992, 1994, 1999; Maoz et al. 1993; Dietrich et al 



ichcrt et al. 1994; Stirpe et al. 1994; AUoin et al 



1991 



NGC 4151 (Cre nshaw et al. |1996|; K aspi et al. |1996[ War- 
wick et al. 1996t Edelson et al. 1996), Fair all9 (Rodriguez- 



1993 



Korista et al.|1995| ; Chiang et al.|2000D , NGC 378 3 (Re 



1995) 



Pasc ual et al. |1997| ; Santos -Lleo et al. |l997|) , 3C 390.3 (Leighly 
et al. |l997| ; Dietrich et al. |l998|; O'Brien et a l |l998|) , and 
NGC 7469 (Wanders et al. |l997| ; CoUier et al. |l998|) ' 



A broad-line region (BLR) size of the order of less 
than a few light weeks is indicated for Seyfert 1 galax- 
ies by the correlated variations of the broad emission-line 
flux and of the optical/ultraviolet continuum. It is gen- 
erally assumed that the central supermassive black hole 
is surrounded by an accretion disk. Such an accretion 
disk is the probable origin of a significant fraction of the 
broad emission- line flux (e.g. Laor & Netzer 1989; Dumont 
& Colhn-Souffrin |1990t Halpern |1990| Zheng, Veilleux, 



& Grand i |1991| ; Hubeny et al. |2000D . Recently, Kaspi et 
al. ( 2000D published results of a long term monitoring cam- 
paign of quasars. Including the results of Seyfert 1 galaxies 
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(Wandel et al. 199£ ) they suggested a relation between the 
BLR size and the optical luminosity given by r^-^^L'''^. 

The Seyfert 1 galaxy NGC 5548 has been continuously 
monitored in the optical since late 1988 by the interna- 
tional AGN watch consortium (Peterson et al. 1999; cf. 



of the EUV wavelength range (Chiang et al. 200C). We 
also discuss the shape of the Ha and H/3 line profiles. We 
summarize our results in section 4. 



AUoin ( |1994D for an AGN watch overview). In June 1998 
this prominent Seyfert 1 galaxy was targeted for a coordi- 
nated intense monitoring campaign using EUVE, ASCA, 
and RXTE to study the high energy continuum emission 
and its temporal characteristics (Chiang et al. 200C). Chi- 
ang et al.'s observations indicate that the variations at 
~0.2keV (EUV) appear to lead similar variations at ener- 
gies larger than ~1 keV by 3-8 hrs. This was unexpected 
as it was generally assumed that correlated variations of 
the EUV, UV, and optical emission would all be due to 
reprocessed higher energy radiation. 

Since this campaign provided the rare opportunity to ac- 
cess the high energy continuum variations with especially 
high temporal sampling, we organized a simultaneous ground-column (9), the spectrograph sht width (in the dispersion 



2. Observations and Data Analysis 

Spectroscopic and broad-band photometric measurements 
of NGC 5548 were obtained by several groups in June 
1998. Table 1 gives a brief overview of the individual par- 
ticipants who recorded the data we report here. Each 
group (column 1) was assigned an identification code given 
in column (2) which is used throughout this paper. This 
code was based on that used for the ongoi ng NGC 5548 
monitoring campaign (Peterson et al. 1999 ). Column (3) 
gives the aperture of the telescope used. Columns (4) - 
(8) list the sizes of the focal-plane apertures used in vari- 
ous observations; fixed instrument apertures were used. In 



based campaign for the optical wavelength domain. 

Assuming part of the BLR fiux is emitted from the 
accretion disk, variations on timescales of one day or even 
less should be detected (Stella 1990 ). For this model of a 
relativistic Keplerian disk it was shown that weak narrow 
features should drift across the emission line profile. These 
weak structures are expected to start at high velocities 
i.e. in the outer profile wings, and move towards the line 
center to be due to the longer time delay at larger radii of 
the disk in response to the variable continuum emission. 
The timescale of the shift depends on the mass of the 
central black hole and it is for NGC 5548 of the order 
of less than several days (cf. Stella 1990). The detection 
of such structures would be a strong indication for the 
presence of an accretion disk in the innermost region of 
an AGN. 

Generally, current monitoring campaigns have not pro- 
vided the necessary temporal and spectral resolution for 
detecting such substructures in the broad emission line 
profiles. However short term monitoring campaigns had 
been undertaken for NGC 4151 (Xanthopoulos & DeR- 
obertis 11991 



wick et al. 



Crenshaw et al. 1996; Kaspi et al. 1996 



1996 



War- 



Edelson et al. 1996) and for a small sample 
of AGN (e.g. NGC 5548, NG C415 1, 3C 390.3, Arpl02B, 
Mkn6; Eracleous & Halpern 1993). As yet no significant 
short timescale broad emission line flux variation or profile 
variability has been detected. 

In this paper, we present the results of the optical 
photometric and spectroscopic observations that were ob- 
tained in June 1998 simultaneous with the high energy 
campaign (Chiang et al. 200C). In section 2 we describe 
the optical observations and outline intercalibration pro- 
cedures by which a homogeneous set of photometric and 
spectroscopic measurements can be achieved. In section 3 
we present measurements of the broad-band flux as well 
as of the broad Balmer emission-line flux. We compare 
our results with the results of the simultaneous campaign 



direction) and extraction width (cross-dispersion dimen- 
sion) of the spectra are listed. 

A complete log of the photometric and spectroscopic 
observations is given in Table lA and 2 A (available in 
electronic form at the CDS). 



2.1. Optical Photometry 

Photometric observations of NGC 5548 were obtained by 
Gaskell et al. (sample Q) and Merkulova (sample D). The 
brightness of NGC 5548 was determined with respect to 
stars in the held of NGC 5548 using the photometric se- 
quence (s tars 1 and 2) deflned by Pension, Pension, & 
Sandage ( 1971 ), and star CI and C defined by Lyuty 
( 1972| ). The bright star «1' east to the galaxy located at 
P. A. = 240° at a distance of «9' with respect to NGC 5548 
also was used. This star is referred in the HST guidestar 
catalogue as gsc 0201001062 (RA = l?'* 17™ 29?66, Dec = 
-|-25°03'10'.'5 (2000.0)). Star CI in Lyuty is identical to 
star 4 in Pension, Pension, & Sandage. The U, B, V magni- 
tudes of CI and 4 are identical within the errors and differ 
by only ±0.02 mag (B,V) and ±0.17 mag (U) (cf. Tab. 2). 
Star C is the bright star close to star 3 in Pension, Pen- 
sion, & Sandage located approximately lf3 southeast of 
it. 

The photometric observations of NGC 5548 of sample 
D obtained with the 1.25 m telescope of the Crimean As- 
trophysical Observatory were recorded using a five channel 
version of the Double Image Chopping Photometer - Po- 
larimeter (Piirola 1973). The observations were made si- 
multaneously in five colors using dichroic filters to split the 
light into five spectral bands. The resulting passbands are 
close to the standard Johnson UBVRI photometric sys- 
tem with effective wavelengths at 3600 A, 4400 A, 5300 A, 
6900 A, and 8300 A, respectively. A diaphragm with two 
equal apertures (15" in diameter) was used in the focal 
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Table 1. Overview of observations; the aperture size is given in units of arcsec 



Source 


Code 


Tel. 


Photometry Aperture 


Spectroscopy 






[m] 


U 


B V R I 


Aperture 


Mt. Hopkins Observatory 


C 


1.5 






3x4.6 


Crimean Astrophysical Obs. 


D 


1.25 


15 


15 15 15 15 


— 


Calar Alto Observatory 


G 


2.2 






2.06x2.61 


Special Astrophysical Obs. 


LI 


1.0 






8x6.6,8x19.8 


Special Astrophysical Obs. 


L2 


6.0 






2x6 


Univ. of Nebraska 0.4 m 


Q 


0.4 




— 8 — — 




Shajn reflector, Crimean Obs. 


W 


2.6 






3x11 



Table 2. Adopted magnitudes for standard stars 



band 


1" 


2" 


4,01'' 


c 





U 


14.52 


16.24 


13.23 


11.17 




B 


14.45 


16.07 


12.78 


10.96 




V 


13.80 


15.38 


11.90 


10.47 


11.29 


R 






11.17 


9.96 




I 






10.64 


9.57 





"from Pension, Pension, & Sandage (1971), estimated uncer- 
tainty 0.02 mag. 



''from Lyuty ( 1972 ), estimated uncertainty 0.01 mag. ( U) of CI 
0.02 mag. 



plane of the 1.25 ni telescope. The distance between the 
aperture centers was 26". A rotating chopper alternately 
closes one of the apertures, leaving the other free, thus 
the photocathode is illuminated alternately by the galaxy 
(or the star) and the sky apertures. Because the centers 
of the apertures are closely spaced, background observa- 
tions were also obtained at a distance of about 7' from the 
galaxy nucleus to correct the observational background 
for a contribution of the outer regions of the galaxy. The 
telescope is fully automated and an autoguider was used. 
Observations of the nucleus of NGC 5548 were made on 
nights with good atmospheric conditions — when the es- 
timated seeing was in the range of 1" to 3". The positional 
accuracy by the autoguider during observations is better 
than 20% of the estimated seeing, i.e. better than 0'.'6 in 
the worst case. The galaxy nucleus was also positioned in 
the aperture using the autoguider, with similar position- 
ing errors. The conventional technique of differential mea- 
surements was applied. Two comparison stars (label led as 
CI and C) were taken from the list of Lyuty ( 1972 ). Ob- 
servations were performed against a comparison star CI; 
and the second star, C, was used to check the results of 
the first comparison. In addition to stars CI and C, sec- 
ondary UBVRI standards by Neckel & Chini ( |198C| ) were 
used for an absolute calibration. The measurements were 
carried out by observing in the following sequence, C-Ci- 
sky-AGN-sky~Ci~C-sky. The time resolution was about 
3.5 minutes. During a single observation, 8 integrations 
of 10 sec each were made. Photon statistics (corrected for 
sky background) were applied to calculate photometric 



errors, which generally were the same as the rms errors 
obtained by averaging the 8 integrations. Generally, the 
atmospheric seeing during the observations obtained with 
the 1.25m CAO telescope averaged 2". It was better for 
two nights (June 27/28 and June 30/July 1) and worse for 
the night June 22/23 (2-3 "). 

The V-band observations of Gaskell et al. (sample Q) 
were obtained with the University of Nebraska 0.4 m tele- 
scope in Lincoln, Nebraska. The frames were recorded with 
a ST-7 CCD camera. The flux of NGC 5548 was deter- 
mined for an 8" diameter aperture. The sky was mea- 
sured in an annulus of radius 15 to 20 ". The V-band 
magnitude was derived relative to the comparison stars 1 
and gsc0201001062 i n the field of NGC 5548 (cf. Pension, 
Pension, & Sandage 1971). The exposure time of an in- 
dividual frame was 3 minutes. Generally, to increase the 
signal-to-noise ratio of the flux measurements, the V-band 
magnitudes of three subsequently recorded frames were 
averaged. The seeing during the observations was typi- 
cally 4". There was no correlation detected between the 
derived magnitude of NGC 5548 and the seeing. 

To combine the broad band fiux measurements of sam- 
ple D and Q, the apparent magnitudes were transformed 
into fluxes. The conversion was performed using the fol- 
lowing equations (Allen 1973| ; Wamsteker 1981): 

logFA(3600A) = -0.4771,7-8.361, (1) 
log Fa (4400 A) 
log Fa (5500 A) 



log Fa (7000 A) 
log Fa (9000 A) 



-0. 4771^7 
-0.4777b 
-0.4777y 
-0.4777fl 
-0.4777/ - 



-8.361, 

- 8.180, 
-8.439, 

- 8.759, 
9.080. 



The broad band V-fluxes of sample Q were scaled to 
the V-band fluxes observed for sample D. The additive 
scaling factor was determind by comparison of observa- 
tions which were separated by less than 2 hours. We de- 
rived a scaling factor of (2.944±0.524) 10~^^ erg s^^ cm"^ A^^ 
based on 5 epochs which corresponds to Amag = 0.20±0.02. 
The difference is caused by the significantly different sizes 
of the apertures used to measure the fiux of NGC 5548 
(Tab. 1). The resulting U,B,V,R,I broad-band light curves 
are given in Table lA (available in electronic form at the 
CDS) in units of 10^^^ ergs^^ cm~^ A~^ and are displayed 
in Figure 1. 
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2.2. Optical Spectroscopy 



1.1 


1 1 1 1 1 1 

; V-band 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1.0 


- star 0,C 


0.9 


_l 1 1 1 1 L. 





T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

I-band 



I ' 1*1' 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — h 
R-band 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — h 



V-band 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — h 



B-band 

I 



H — I — I — h 



H — I — I — h 



H — I — I — I- 



H — I — I — I- 



U-band 



H — h 



J I I I I I I I I I I I I I I I I I I I I L 

970 980 990 1000 

JD - 2450000 

Fig. 1. Optical broad-band light curves for NGC 5548. 
Fluxes in broad-band U , B,V , R, and / are in units of 
10~^^ erg s~^ cni~^ A~^. In the top panel the normalized 
V-band flux of star (-1-) and C (x) is shown which was 
contant within ^^0.57% . 



The flux calibration of AGN spectra can be accomplished 
in several ways. In variability studies, it has become com- 
mon practice to normalize the flux scale to the fluxes 
of strong forbidden narrow emission lines which are as- 
sumed to be constant over timescales of at least several 



decades (e.g. Peterson 1993 ). This assumption is justified 
by the large spatial extent and low gas density of the 
narrow- line region (NLR), but is only valid if the aper- 
ture used is larger than the NLR. Light travel-time effects 
and the long recombination timescale (tj-oc ~ 100 years for 
rte « 10'^ cm^'^) damp out short timescale variability. 

As a major sources of uncertainty in inhomogenous 
samples of variability data is the use of different instru- 
mental settings, it is important to t ake aperture effects 
into account (Peterson & Collins 1983 ). The seeing-dependent 
uncertainties which are introduced by the aperture geom- 
etry can be minimized by using large apertures. It has 
been shown that apertures of 5" x 7'.' 5 can reduce seeing- 
dependent photometric errors to no more than a few p er- 
cent in the case of nearby AGNs (Peterson et al. 1995| ) . 

The spectra of NGC 5548 obtained at the 1 m tele- 
scope of SAO RAS (sample LI Table 1) were extracted 
for a slit aperture of 8'.'0 x 6'.'6 and 8'.'0 x 19'.'8, respec- 
tively. For these large apertures seeing losses can be ne- 
glected. To measure the F([0 III] 5007 A) a linear pseudo- 
continuum was fitted beneath the [O III] emission line and 
the flux was measured for the wavelength range AA5062- 
5108 A. For the extraction windows used for sample LI 
we determined a flux of F([OIII] 5007 A) = (5.07±0.05) 
10-"ergs-icm-2 (8'.'0 x 6'.'6) and F([OIII] 5007 A) = 
(5.91±0.05) 10-i3gj.gg-icm-2 (gz/g x 19'.'8), respectively. 
These values are in good agreement with F([0 III] 5007 A) 
measured for NGC 5548 during the last 10 years (cf. Peterson 
et al. 1999). For large apertures of comparable size (e.g. 
8'.'8 X 12'.'8) F([0 III] 5007 A ) = 5 1 to 5.9 lO^i^ergs-i cm^^ 
is mentioned (Peterson et al. |l992| ). To be constistent with 
the H/3 line flux measurements of NGC 5548 obtained since 



1988 (Peterson et al. [l999| ) we adopted F([OIII] 5007 A) 

-13^ 



(5.58±0.27) 10-"ergs-icm- 



2.3. Intercalihration of the Spectra 

Since the spectra of the various samples were taken with 
different instruments in different configurations, it was 
necessary to intercalibrate them to a common flux level. 
As we have shown above, the [O III] 5007 line flux was 
constant to better than 1% for large extraction windows. 
We can safely use the narrow emission lines such as the 
flux of [O III] 5007 as standards. In order to avoid any 
wavelength-dependent calibration errors, each spectrum 
was scaled in flux locally over a limited wavelength range 
prior to measurement. The H/3 spectral region was scaled 
with respect to the [O III] 4959,5007 line fluxes. The Ha 
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Table 3. Scaling factors for photometric and spectro- 
scopic subsets 



sample 


add. constant 


point source 


extended source 




V-mag." 


correction if^ 


correction G'' 


C 


— 


0.972±0.022 


2.02±0.33 


D 


1.000 






G 




0.865±0.022 


3.36±0.54 


LI 




1.00 


0.00 


L2 




0.945±0.025 


2.60±0.53 


Q 


2.944±0.524 






w 




0.950±0.028 


1.52±0.32 



° relative to sample D 
''relative to sample LI 



region was scaled vifith respect to the flux of the [O I] 6300, 
[Nil] 6548,6583, and [S II] 6716,6731 emission lines. The 
spectra were intercalibrated using the method described 
by van Groningen & Wanders ( 1992| ). This procedure cor- 
rects for different flux scales, small wavelength shifts, and 
different spectral resolutions by minimizing the narrow- 
line residuals in difference spectra formed by subtracting 
a "reference spectrum" from each of the observed spec- 
tra. The rescaled spectra were used to derive integrated 
emission-line fluxes. 



3. Results 

3.1. Light Curves 

The broad band flux variations in U,B,V,R,I are dis- 
played in Fig. 1. In addition to the broad band flux vari- 
ations the light curve of the measured flux of the com- 
parison stars C and (cf. Tab. 2) is shown. The normal- 
ized flux remained constant within ^--^0.57% . The broad 
band flux of NGC 5548 decreases throughout the cam- 
paign in June 1998. This trend is clearly visible in the U- 
and B-band, and the ratio of the flux levels observed in 
early and late June becomes smaller for increasing wave- 
length. In the I-band the flux decay is negligible. The 
merged V-band light curve of sample D and Q shows 
different internal scatter of the measurements recorded 
during individual nights (Figs. 1,2). This difference might 
be caused by different seeing conditions of the observa- 
tions. Merkulova (sample D) measured the brightness of 
NGC 5548 with an 15" aperture under a typical seeing of 
1-3". The seeing amounts to «4" for the measurements 
Gaskell et al. (sample Q) observed for an 8" aperture. 
Modelling the host galaxy of NGC 5548 Romanishin et 
al. ( |1995| ) showed that the light loss of the AGN is of the 
order of 10% for sample Q and only slightly larger for the 
host galaxy. However, the different scatter may be also in- 
trinsic. The temporal resolution of the samples amounts to 
<4 minutes (D) and nearly 16 minutes (Q), respectively. 
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987.35 987.40 987.45 987.50 987.55 

JD - 2450000 

Fig. 2. The strong variation in the broad band fluxes as 
measured on June 22 (JD=2450987). In all 5 broad band 
fluxes a rapid increase and subsequent decay was detected 
within w90 minutes. In the top panel the normalized V- 
band flux of star (+) and C (x) is shown. 



The scatter of the brightness during individual nights is 
~5% for sample D. This amplitude is of the order which is 
expected for microvariability. For several Seyfert 1 galax- 
ies (NGC 3516, NGC 4151, NGC 5548, NGC 7469) detec- 
tion of microvariability was claimed on timescales of less 
than 10 minutes (cf. Lyuty et al. 198£; Dultzin-Hacyan et 
al. 
al. 



1992; Jang & Miller 1997; Noble et al. 1997; Carini et 



1998; Welsh et al. 1998; Ghosh et al. 200C; Merkulova 



pOOOD and amplitudes of the order of 15% to 20% . How- 



ever, it seems that onset and disappearance of microvari- 
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Table 4. Integration limits 



ICdL Lli. C 


ur£n7f> 1 on o'^" h n o'p 

WCL V Lll 1 CLll tic 




[A] 


Hq A6563 


6500-6800 


H/3A4861 


4850-5020 


Fa (5100 A) 


5185-5195 



ability follow a random process since several attempts to 
detect or to confirm rapid optical variations of Seyfert 1 
ga laxies failed (cf. Dultzin-Hacyan et al. 1993 ; Petrucci et 
al. 199£ ). Hence, it might be that rapid variations can not 
be detected on timescales larger than 15 minutes. 

The most obvious structure in the broad band light 
curves is a flare like event. This variation was detected on 
JD = 2450987 (June 22). The observations recorded during 
this night are shown in more detail in Fig. 2. In all 5 broad 
band measurements a short intense increase of the flux 
level occured. The broad band flux increased by 20±2% 
within w 30 minutes and decreased within «60 minutes 
to the flux level of the beginning of this event. The dura- 
tion of this flare lasted only ^90 minutes. This remark- 
able flux increase can not be due to temporal variations 
of a standard star which was used for relative flux cali- 
bration. The flux calibration is based on two stars which 
were in cross checked with secondary standard stars taken 
from Neckel & Chini ( |198QD in addition. As can be seen 
in Fig. 1 and 2 the used comparison stars were constant 
within less than 0.57% . The standard deviations of the 
primary star used for calibration during this night amount 
to 0.018,0.014,0.009,0.023,0.001 mag in U,B,V,R,I. 

Broad emission-line fluxes were integrated over a wave- 
length range of AA4850-5020 A for H/3 and AA6500-6800 
A for Ha. A local linear continuum fit was interpolated 
under each emission line. In the case of the H/3 region, the 
continuum was defined by the flux measured in two win- 
dows (lOA width) at 4845 A and 5170 A in the observed 
frame (H/3) and at 6350 A and 6965 A (Ha). No attempt 
was made to correct any of the measured emission-line 
fluxes for their respective narrow-line contributions. The 
optical continuum flux Fa (5100 A) was determined as the 
average value in the range 5185-5195 A (cf. Tab. 4). 

The measured H/3 and Ha emission line fluxes were ad- 
justed to take into account the different slit sizes (Tab. 1). 
The measurements of both H/3 and Ha were normalized 
by a similar factor 



FHfi = F[o III] 



'^[o ///] 



(2) 



obs 



The continuum fluxes Fa (5100 A) were then adjusted 
for different amounts of host-galaxy contamination (see 
Peterson et al. 1995 for a detailed discussion) by 

"Fa (5100 A) 



Fa (5100 A) ^ipFi 



[OH I] 
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Fig. 3. Light curves for the emission lines Ha, H/3, and the 
optical continuum flux Fa (5100 A). The vertical scale is in 
units of 10""'^'^ erg s~^ cm~^ A~^ for the lines and 10""'^^ erg 
s^^ cm^^ A~^ for the continuum. 



where F^qiii] is the adopted absolute flux of F([OIII] 
5007 A), the quantity in brackets is the observed contin- 
uum to F([OIII] 5007 A) flux ratio measured from the 
spectrum, and G is an aperture-dependent correction for 
the host-galaxy flux. The sample LI which uses a relatively 
large aperture (8'.'0x 6'.'6) was adopted as the standard 
(i.e., (/3=1.0, G = by definition), and other data sets 
were merged progressively by comparing measurements 
based on observations made during the same night. Note 
that this means that any real variability that occurs on 
timescales this short will be somewhat suppressed by the 
process that allows us to combine the different data sets. 
The scaling factor tp and the additive scaling factor G for 
the various samples arc given in Table 3. 

The resulting light curves of the optical continuum 
Fa (5100 A) and of the H/3 and Ha emission line fluxes are 
displayed in Fig. 3 and presented in Table 3A (available 
in electronic form at the CDS). The optical continuum 
light curve is very similar to the broad-band V light curve 
(cf. Figs. 1,3). The difference in the ffux levels is probably 
due to the larger aperture used for the broad band flux 
determination relative to that used for the optical contin- 
uum fluxes. 
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feature a^^^^ a^^^ feature a^^^^ a^^^ feature N interval interval" 

U 0.008 0.021 I 0.006 0.024 I'^^y^l t'^^y^l 



B 0.008 0.020 Fa (5100 A) 0.020 0.016 


u 


120 


0.24 ± 1.53 


0.03 ± 0.15 


V 0.010 0.020 H/3 0.022 0.010 


B 


120 


0.24 ± 1.53 


0.03 ±0.15 


R 0.004 0.016 Ha 0.040 0.025 


V 


178 


0.20 ±0.87 


0.15 ±0.47 




R 


120 


0.24 ± 1.53 


0.03 ±0.15 


"observational uncertainty based on uncertainties assigned to 


I 


120 


0.24 ± 1.53 


0.03 ±0.15 


individual points. 


Fa(5100A) 


52 


0.59 ±1.96 


0.32 ± 0.36 


'mean fractional uncertainty based on point-to-point differ- 


H/3 


52 


0.59 ±1.96 


0.32 ± 0.36 


ences between closely spaced (i.e., At < 0.05 d) measurements. 


Ha 


38 


0.81 ±2.28 


0.45 ± 0.43 



A final check of the uncertainty estimates was per- 
formed by examining the ratios of all pairs of photometric 
and spectroscopic observations which were separated by 
0.05 days or less. 

There are more than a dozen independent pairs of mea- 
surements separated by less than 0.05 d. For Ha, the un- 
certainty estimate is dominated by the spectra taken at 
Calar Alto Observatory, while for H/3 more than 20 pairs 
could be used. The dispersion about the mean (unity), 
divided by \/2, provides an estimate of the typical uncer- 
tainty in a single measurement (aest). The observational 
uncertainties (cobs) assigned to the spectral line flux mea- 
surements were estimated from the error spectra which 
were calculated by the intercalibration routine, as well as 
from the signal-to-noise ratio within the spectral range 
near the individual emission lines. For example the mean 
fractional flux error is (Tobs=0.040 for the Ha line, and the 
average fractional uncertainty from the internal statisti- 
cal estimate is (Test = 0.025. For the H/3 line the values 
are tTobs=0.022 and (Test ~ 0.010 which implies that the 
error estimates for both lines are probably quite good. 
Generally, the estimated errors ((Tost) are of the same or- 
der as, but slightly smaller than, the observational un- 
certainties ((Tobs) derived directly from multiple measure- 
ments (Tab. 5). The large difference of (Tobs and (Tost for 
the broad band photometric fluxes is mainly due to the 
measurements obtained on JD=2450987 (June 22), which 
we expect is due to a real variation. 

The average interval between measurements of the com- 
bined broad band light curves is about 0.24 ± 1.53 days 
for U, B, R, I, and for V 0.20 ± 0.87 days. However, 
if the large gaps (>5days) are not taken into account, 
the sampling rate drops to 0.03 ± 0.15 days {U,B,RJ) 
and 0.15 ± 0.47 days (V) (Tab. 6). The sampling of the 
Ha, H/3, and the optical continuum flux F^ (5100 A) can 
be obviouly devided into two regimes. Taking into ac- 
count all measurements the sampling is about Ai(H/?) 
= 0.59±1.96days for H/3 and Fa(5100A). Neglecting the 
large temporal gap of nearly two weeks at the beginning 
of the monitoring campaign (full moon). At is Ai(H/?) = 
0.32±0.36days. For the Ha light curve the correspond- 
ing sampling rates are Ai(Ha) = 0.81 ± 2.28 days and 
Ai(Ha) = 0.45 ± 0.43 days, respectively. 



^neglecting temporal gaps larger than 5 days 



The broad Ha and H/3 emission lines, the optical con- 
tinuum flux Fa (5100 A) and the U,B,V,R,I broad band 
fluxes appear to exhibit small amplitude variations on 
timescales of days (Figs. 1-3). The broad band light curves 
U ^By ^ as well as the F^ (5100 A) flux show a decreasing 
flux from the beginning of this campaign to the end. The 
broad band flux at longer wavelengths [R and /) appeared 
to be nearly constant. 

The Ha and H/3 light curves show weak evidence for 
small amplitude variations on timescales of only a few 
days. However, this is statistically not significant; within 
the errors the Ha and H/3 emission line flux can be re- 
garded as nearly constant. This can be seen by deriving 
the variability parameter F„ar. 

The variability parameters i*Var and i?max have been cal- 
culated for the broad-band flux variations as well as for the 
broad emission lines and the optical continuum (cf. Clavel 
et al. 1991; Rodriguez-Pascual et al. 1997). The quantity 
^max is the ratio of the maximum to the minimum flux. 
The quantity i^var is an estimation of the fluctuations of 
the intrinsic variations relative to the mean flux. There- 
fore, the rms of the light curves has been corrected with 
respect to the uncertainties introduced by the observa- 
tions. Values of i?max and i^var for the broad-band mea- 
surements are given in Table 7. The value of Fvar for the 
optical contimmm Fa (5100 A) agrees well with Fvar for the 
V-band. Furthermore, for the broad band variations and 
the optical continuum we reproduce the well known trend 
that the amplitude of the variation decreases with increas- 
ing wavelength (Tab. 7). The broad band fluxes and the 
Fa (5100 A) continuum flux were not corrected for the host 
galaxy contribution since this fraction is only of the order 
of less than 10% (Romanishin et al. 1995). The values 
of ¥yar determined for the Ha and H/3 emission line flux 
variations indicated that the Balmer line flux can be taken 
as constant within the errors, even if one is tempted to 
glimpse hints of small amplitude variations on timescales 
of a few days (cf. Tab. 7, Fig. 3). 

For comparison we also calculated Fvar and Rmaa; for 



the rapid variations detected by Chiang et al. (2000) in 
the extrem ultraviolet domain (Tab. 7). In Fig. 4 we show 
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Table 7. Variability statistics of the light curves; con- 
tinuum flux and broad band fluxes in units of 10~^^erg 
s~^cm^^A^^; line fluxes in units of 10^^^ erg s^^cm"^. 
The EUV flux is given in units of counts per second. 



feature 


mean flux 


rms flux 




J- var 


U 


24.80 


2.37 


1.48 


0.0865 


B 


20.08 


1.56 


1.44 


0.0697 


V 


17.76 


0.91 


1.38 


0.0421 


R 


17.54 


0.78 


1.34 


0.0406 


I 


13.15 


0.69 


1.44 


0.0557 


Fa (5100 A) 


14.4 


0.9 


1.36 


0.040 


H/? 


10.5 


0.2 


1.11 


-0.005 


Hq 


38.7 


1.2 


1.13 


-0.006 


EUV 


0.0993 


0.0376 


59 


0.2035 



1.2 



1.0 



o.g 



E 2.5 
° 2.0 
1.5 
1.0 
0.5 
0.0 



— 1 1 1 1 1 1 1 1 1 r- 
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H 1 1 h 
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Fig. 4. The normalized EUVE light curve from the moni- 
toring program for NGC 5548 (Chiang et al. 2000) is shown 
in the bottom panel. In the top panel the normalized vari- 
ations of Fa (5100 A) are displayed. 



the normalized light curve of Fa(5100 A) and of the EUV 
variations (kindly provided by Chiang and collaborators). 
The Fa (5100 A) continuum flux show the same overall 
trend as the EUV variations. But in the optical contin- 
uum the rapid variations visible in the EUV are smeared 
out or the amplitude of the variations at A=5100A is so 
small that it is hidden by the measurement uncertainti- 
ties. The disappearance of small scale variations with in- 
creasing wavelength was al so ob served for NGC 4151 in 
Dec. 1993 (cf. Edelson et al.|l996D. 



3.2. Time-Series Analysis 



The emission lines are expected to change in response to 
variations in the far-UV continuum, primarily in response 
to variations at the unobservable wavelengths just short- 
ward of 912 A. Generally, it is necessary to assume that 
the observable continuum can approximate the behavior 
of the ionizing continuum and that the shortest observed 
UV wavelengths provide the best observable approxima- 
tion to the "driving" (ionizing) continuum. However, for 
the first time we can make use of a well-sampled simul- 
taneous X-ray light curve on short timescales (Chiang et 
al. 200C ) . The EUVE light curve and the variations of the 
optical continuum show very similar variability pattern 
(Figs. 3 and 4). 

However, the Ha and H/3 emission line fluxes show only 
marginal indications for variations. Within the errors the 
emission line flux was constant within 5% which is also 
indicated by Y^ar (cf . Table 7) . In spite of this we calcu- 
lated cross-correlation functions. We used the extrem ul- 
traviolet light curve measured with EUVE for the driving 
continuum kindly provided by Chiang and colla borat ors. 
We computed the ICCF (cf. White & Peterson |l994D i n 
the formalism of the so called "local CCF" (Welsh |l 9991) . 
This approach takes into account the influence of the bias 
t o und erestimate the delay as has been shown by Welsh 
( 1999| ) . To correct for the influence of low frequency power 
in the flux variations we calculated a linear trend which 
was subtracted. The removal of a linear trend helps to re- 
duce the bias of the observed delay towards values smaller 
than the real delay (Welsh 199E ). We also applied the 
cross-correlation analysis to the detrended light curves, 
and we calculated the discrete correlation function (DCF) 
for comparison. Uncertainties in the ICCF results for the 
cross-correlation maxima and centroids were computed 
through Monte Carlo techniques (cf. Peterson et al. 1998). 

The sampling characteristics of each of the light curves 
are given in Table 8. The name of the feature and the total 
number of points, iV, in the light curves are given that 
were used in computing the auto- and cross-correlation 
functions. The width (FWHM) of the AGE based on the 
original and of the detrended light curves are given as well. 
The corresponding widths of ICCF computed by cross- 
correlation with the EUV continuum at ~ 0.2 keV is listed 
too. 

The results of the cross-correlation analysis are given 
in Table 9. The first column indicates the "responding" 
light curve (i.e., the light curve that is assumed to be 
responding to the driving light curve). The second and 
third column provide the peak value of the correlation 
coefficient r^ax for the ICCF (original and detrended). 
The position of the peak of the cross-correlation functions 
''"max was measured by determing the location of the peak 
value of the ICCF and DCF; these values are given in 
the next four columns, respectively. The centroids Tccnt of 
the ICCF of the original and detrended light curves were 
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Table 8. Sampling characteristics and full width at half 
maximum (FWHM) of ACF and CCF features. The 
FWHM is given in units of days 



delay T[d] 



Fig. 5. Cross-correlation functions for U, B, V, R, I fluxes, 
Fa (5100 A), and the H/3 and Ha emission line (from bot- 
tom to top). The ICCF computed with the original light 
curves (dashed line) and the ICCF based on the detrended 
light curves (solid line) are displayed. The average ICCFs 
which were obtained by estimating the uncertainties of 
the delays are shown as short solid line but restricted to - 
10< T <10d. The discrete points show the DCF computed 
for each detrended light curve. 



feature 


N 


FWHM(ACF) 


FWHM(ICCF) 






ong. 




orig. 




U 


120 


6.8 


3.3 


7.1 


5.0 


B 


120 


5.7 


4.2 


6.6 


5.9 


V 


178 


1.7 


0.2 


9.3 


2.2 


R 


120 


1.8 


1.6 


4.1 


5.9 


I 


120 


0.2 


0.2 


3.6 


3.7 


Fa(5100A) 


52 


30 


0.4 


6.1 


3.3 


H/3 


52 


1.7 


0.6 


4.0 


2.4 


Ha 


38 


1.2 


0.2 


1.1 


1.4 



computed using the points in the cross-correlation func- 
tion with values greater than 0.6r,„ax- The error estimate 
for the position of the cross-correlation peak Ar^ax and 
the ICCF centroid Arcont for the original and detrended 
light curves are also given (cf. Tab. 9 columns l^Torig. and 

Ardetr.). 

The broad-band [U , B, V, R, I) and the optical contin- 
uum flux variations appear to be simultaneous within the 
errors, relative to the EUV continuum variations at ~ 0.2 
keV (Fig. 5). Taking into account the small variability am- 
plitude expressed by Fvar (Tab. 7), the temporal coverage 
of the time series of 30 days only, with dense temporal 
sampling solely for the second half, the small CCF ampli- 
tudes, and the uncertainty of the location of the peak and 
centroid of the CCFs of «3 days (Tab. 9) no reliable delay 
was detected for the broad emission lines of Ha and H/3 
in response to the EUV variations (Fig. 5). 



3 3.3. Mean and Root-Mean-Square Spectra 



We calculated mean and root-mean-square (rms) spectra 
from the flux-scaled spectra for the Ha and H/3 regions. 
The rms spectrum is useful for isolating variable parts of 
line profiles. 

The mean and rms spectra of the samples C, LI, L2, W are 
identical within less than 3% (H/3) and 1.5% (Ha). The 
mean and rms spectra of the Ha and H/3 region are pre- 
sented for the spectra of subsample C (Fig. 6). It provides 
19 epochs during the studied period and a daily sampling 
with homogeneous settings for the second half of June 
1998. The rms spectrum of H/3 and Ha indicates contin- 
uum variations of the order of less than 5% . Furthermore, 
weak broad emission line flux variations are superimposed 
but the amplitude is less than 2% . Thus, these features 
in the rms-spectra again should be taken only as weak in- 
dications for emission line flux variations on timescales of 
days. However, the width of the broad Hell 4686 emission 
line in the rms spectrum is similar to that detected by 



Peterson et al. (2000) for the narrow-line Seyfert 1 galaxy 
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Fig. 6. The mean and root-mean-square (rms) spectra the of the H/3 (top) and Ha (bottom) line region are shown 
for sample C. The vertical scales are arbitrary. The rms spectrum of H/3 and Ha provide some indication for emission 
line variations but only of the order of less than 2% . Remarkably the broad feature at the location of Hell 4686 in 
the rms spectrum of the H/3 wavelength region. It is significantly broader than those of the H/? or H7 line. 

Table 9. Results of the cross-correlation analysis for the orginal and detrended light curves. 





r 


max 






[days] 




Tcent 


[days] 


error, est. 


feature 


orig. 


detr. 


ICCF 


DCF 


ICCF 




days] 








orig. 


detr. 


orig. 


detr. 


orig. 


detr. 




U 


0.86 


0.57 


0.0 


0.9 


0.0 


1.0 


-0.1 


-0.1 


0.6 


1.1 


B 


0.83 


0.65 


-0.1 


-1.0 


0.0 


1.0 


-0.4 


-0.6 


0.6 


0.6 


V 


0.71 


0.40 


1.0 


1.0 


-2.0 


1.0 


1.4 


1.0 


2.1 


2.1 


R 


0.60 


0.62 


-1.2 


-1.3 


1.0 


1.0 


-0.0 


-0.1 


0.9 


0.8 


I 


0.52 


0.58 


0.5 


0.5 


1.0 


1.0 


0.3 


0.6 


5.2 


2.0 


Fa (5100 A) 


0.83 


0.53 


-0.1 


-0.1 


0.0 


0.0 


0.5 


0.1 


1.4 


3.2 


H/3 


0.56 


0.04 


6.1 


6.5 


6.0 


7.0 


6.1 


7.8 


3.0 


4.5 


Hq 


0.26 


0.29 


9.8 


9.7 


10. 


10. 


9.8 


9.7 


6.0 


2.7 



NGC 4051. The contribution of FeH emission (mulitplets 
37,38) should be negligible since there is no indication of 

variable FeH emission (multiplets 48,49) in the rms spec- 
trum in the wavelength range AA 5200-5500 A. 

The spectra obtained at Calar Alto Observatory have 
a significant higher spectral resolution (R~8kms~^) than 
the spectra of the other samples (R~500kms~^). Hence, 
we used the spectra of sample G to search for small fea- 



tures as suggested by Stella (199C). The H/3 profile was 



filtered with a narrow gaussian curve (FWHM=0.lA) to 
improve the signal-to-noise ratio. The average spectra of 
each night of the 15 epochs obtained at Calar Alto Ob- 
servatory are presented in Fig. 7 as well as the mean spec- 
trum of the samples C,L1,L2, W. The shape of the H/3 
profile is smooth and no narrow features have been de- 
tected during this short period. There is no obvious nar- 
row moving feature visible in the H/3 profile. However, at 
A ~ 4962A there might be a marginal indication of a weak 
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Fig. 7. The average echelle spectra of the H/3 emission hne 
profile as obtained at Calar Alto Observatory are shown 
for each night. For comparison the mean spectrum ob- 
tained from the samples C,L,W with lower spectal reso- 
lution is displayed at the bottom of the panel. The vertical 
dotted lines indicate the location of the marginal feature 
at A ~4962 A as well as the corresponding blueshifted lo- 
cation. 



structure which is visible in the mean spectrum and in 
most of the individual spectra of sample G. It is redshifted 
by --llOOkms"^ and the width amounts to FWHM~6A 
corresponding Av ~350kms^^. We used also the spec- 
tra of the samples with lower spectral resolution to search 
for this weak feature. It is also detectable in the mean 
spectrum based on these samples as can be seen in Fig. 7. 
Unfortunately, the redshift of v=1100kms^^ places this 
weak feature in the red wing of the prominent [N II] 6583 
emission line. Hence, it is not possible to detect this struc- 
ture in the Ha line profile to provide further evidence for 
the existence of this weak feature. 



4. Summary 

The results of a short optical monitoring campaign (June 
1998) on the Seyfert 1 galaxy NGC 5548 are presented in 
this paper. The principal findings are as follows: 

1. The broad-band {U,B,V) fluxes, and the optical con- 
tinuum measured from spectrophotometry iOv(5100 A), 
showed a monotonic decrease of approximately 35 % 
([/), 25% (B), 15% (V), and 30% (Fa (5100 A)) from 
the beginning to the end of the campaign (~ 30 days). 

2. The broad-band (i?,/) fluxes, and the integrated emis- 
sion-line fluxes of Ha and H/3 showed only marginal 
indications for variations on timescales of the order of 
a week with amplitudes less than 5% . 

3. In all broad band fluxes a short flare like event was 
detected during the night of June 22 (JD 2450987) 
which lasted ss90 minutes. 

4. The shape of the variability pattern of the optical con- 
tinuum is very similar to the detected variations in the 
EUV energy range 0.2 keV) with EUVE by Chiang 
et al. ( ^OOOD . But F^ar and R,„ax of the EUV varia- 
tions are significantly larger than Fyar and Kmax of 
Fa(5100A). 

5. The parameter Fvar, which is essentially the rms varia- 
tion about the mean, increased with decreasing wave- 
length for the broad-band measurements. The optical 
continuum Fa (5100 A) and the V- and R-band exhibits 
comparable Fvar- 

The Balmer emission lines Ha and H/3 underwent no 
significant variations as measured by Fvar during this 
campaign within the errors. 

6. The variations of the broad-band and optical contin- 
uum fiuxes are simultaneous with the ~ 0.2 keV vari- 
ations. In spite of the Balmer lines Ha and H/3 showing 
no significant variability, we calculated cross-correlation 
functions. No reliable delay was detected for Ha and 
H/3 with respect to the EUV continuum within the 
uncertainties. 

7. The H/3 line profile was examined for narrow features. 
There is marginal evidence for a weak feature at A ~4962 

A with FWHM~6 A which is redshifted by Av~1100 km s"^ 
with respect to H/3n 



-'narrow • 
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